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Abstract 

Using computer simulations, we show that metastable states still 
occur in two-lane traffic models with slow to start rules. However, 
these metastable states no longer exist in systems where aggressive 
drivers (which do not look back before changing lanes) are present. In- 
deed, the presence of only one aggressive driver in the circuit, triggers 
the breakdown of the high flow states. In these systems, the steady 
state is unique and its relaxation dynamics should depend on the lane 
changing probability p c h and the number of aggressive drivers present 
in the circuit. It is found also that the relaxation time r diverges as 
the form of a power-law : r oc p~h , 0=1- 

PACS. 02.50.-Ey Stochastic processes 05. 45. -a Nonlinear dynamics 
and nonlinear dynamic systems 45.70.Vn Granular models of com- 
plex systems; traffic flow 89.40. +k Transportation 



*e-mail: najemmoussa@yahoo.fr 



1 



Recently, cellular automata (CA) traffic models are used enormously in 
order to understand the complex dynamic behavior of the traffic in roadways 
(see the review pQ). In CA, time and space are discrete. The space is rep- 
resented as a uniform lattice of cells with finite number of states, subject to 
a uniform set of rules, which drives the behavior of the system. These rules 
compute the state of a particular cell as a function of its previous state and 
the state of the neighboring cells. The most popular CA model for traffic 
flow on one-lane roadway is the NaSch model [2]. Despite its simplicity, the 
model is capable of capturing some essential features observed in realistic 
traffic like density waves or spontaneous formation of traffic jams. To de- 
scribe more complex situations such as multi-lane traffic, extensions of the 
NaSch model have been proposed where additional rules are added for lane 
changing cars. 

Barlovic et al [3] found metastable states in their velocity dependent random- 
ization (VDR) which is an extension of the NaSch model. The one-lane VDR 
model belongs to the class of CA models with "slow-to-start" rules. These 
models show an hysteresis effect which is a consequence of the non-unique 
dependence of the flow on the density. The above characteristic behavior 
of traffic flow is also observed in two-lane traffic models. Indeed, Awazu [I] 
showed the appearance of several branches and hysteresis in the relation be- 
tween traffic flow and car density. 

To establish the existence of metastable states, Barlovic et al started their 
simulations of the VDR model from two different initial configurations, the 
megajam and the homogeneous state. The megajam consists of one large 
compact cluster of standing cars. In the homogeneous state, cars are dis- 
tributed periodically with equal constant gap between successive cars (with 
one lager gap for incommensurate densities). If the initial configuration is 
homogeneous, one obtains the upper branch, for some interval of densities 
[pi,P2], in which each car moves freely with maximal velocity (see figure 1). 
This upper branch is metastable with an extremely long life-time. If the ini- 
tial configuration is megajam, one obtains the lower branch which is phase 
separated. The phase separated state consists of a large jam (jammed re- 
gion) and a free-flow region where each car moves freely. It is known that 
the lifetime of the metastable states depends on the system length L. Yet, 
the simulation results indicate that Ap = p\ — p 2 decreases with larger sys- 
tem sizes and is expected to vanish for L — > oo, i.e. the jammed branch is 
stable in that limit. Therefore the non-unique behaviour of the fundamental 
diagram is only observable if finite system sizes are considered (see pQ page 
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93). 

As vehicular traffic usually evolved in multi-lane roads, some interesting ques- 
tion is not yet studied. Does multi-lane version of the VDR model always 
exhibits metastable states? 

The NaSch model with VDR rule is a one-dimensional probabilistic CA which 
consists of N cars moving on a one-dimensional lattice of L cells with peri- 
odic boundary conditions (the number of vehicles is conserved). Each cell is 
either empty or occupied by just one vehicle with velocity v = 1, 2, v max . 
We denote by Xk and Vk the position and the velocity of the kth car at time 
t respectively. The number of empty cells in front of the kth car is denoted 
by dk = Xk+i — Xk — 1 and called hereafter as the gap. Space and time are 
discrete. At each discrete time-step t — ► t+1 the system update is performed 
in parallel for all cars according to the following four subrules : 
R\. VDR, p(vk) = Po for v k — and p{vk) = p for Vk > 0. 
R2. Acceleration, Vk <— min (vk + 1, f max ) . 
R 3 : Slowing down, t> fc <— min (v^, dk)- 

R4. Randomization, vj. <— max (vk — 1, 0) with probability p(vk)- 

R 5 : Motion, the car is moved forward according to its new velocity, Xf. <— 

x k + v k - 

In two-lane traffic models, lane changing of vehicles are performed according 
to some additional rules [5-9]. In this paper, we shall adopt the symmetric 
exchange rules which are defined by the following criteria [B|: 

1. min{v k + l,v max ) > d k 

2. dk,other > dk and d^back > hack 

3. p ch > rand() 

Here d k)0 ther (resp. d^back) denotes the gap on the target lane in front of (resp. 
behind) the car that wants to change lanes. Two different formulas are as- 
signed to the parameter hack- For aggressive drivers, i.e., vehicles which do 
not look back before changing lanes, we choose hack = 0. For careful drivers, 
i.e., vehicles which respect the safety criterion, we choose hack = where 
v b Q is the velocity of the following car in the target lane. Finally, p c h is the 
lane-changing probability and rand() stands for a random number between 
and 1. Hereafter, we shall denote by VDRMjv a the two- lane VDR traffic 
models where iV a represents the number of aggressive drivers present in the 
circuit. 
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The update in the two-lane model is divided into two sub-steps: in one sub- 
step, cars may change lanes in parallel following the above lane changing 
rules and in the other sub-step each car may move effectively by the forward 
movement rules as in the single-lane traffic. 

We performed computer simulations of the two-lane model with the following 
parameters, (p = 0.01, p = 0.7 and v max = 5). The size of the lattice is 
given by L = 1000. Starting from an initial configuration (homogeneous or 
megajam) the system evolved in time steps with respect to the above dynam- 
ical rules. For each simulation run, we discarded some number (tdis) of time 
steps and we performed averages of the flow over t av = 50000 time steps. 
The duration of each run is n t dis + t av " . The procedure is then repeated for a 
number 100 of different realizations of the homogeneous (or megajam) initial 
configurations. The average over all the different realizations gives a mean 
value of the flow. 

Figure 1 illustrated the variation of the flow J, in two-lane VDR traffic mod- 
els, as a function of the density of cars and for different values of the discarded 
time tdis- We noticed that the flow in both lanes are equal since symmetric 
lane changing are considered. First, we shall consider the case where only 
one aggressive driver is present in the circuit. So, if the homogeneous initial 
state is used, a higher branch of the flow is observed for some interval of 
densities [pi, p?\ whenever i dis is small enough. When increasing enough t dis , 
the high branch interval diminished and disappeared completely at certain 
limit of tdis- Notice that this phenomena occurred also for the NS model with 
very small randomization p. In contrast, the upper branch in the fundamen- 
tal diagram of the VDRM does not change when one increases the time tdis 
(Fig. 1). This shows clearly that the hysteresis exist in the fundamental 
diagram of the VDRM model but not in the one of the VDRMi. 
To clarify more the above results, we shall consider the time evolution of the 
flow for some fixed density p = 0.12 (pi < p < p 2 ) and for the homogeneous 
and megajam initial states (Fig. 2). It is shown that in contrast to the 
VDRM , where the homogeneous state is metastable with an extremely long 
life-time, this state does not exist in VDRMi . Yet, in this later, the flow 
corresponding to the homogeneous initial configuration decreases with time 
until reaching the value corresponding to the megajam initial configuration. 
The breakdown of the homogeneous structure in the two lanes is due to the 
occurrence of stopped cars provoked by the abrupt lane changing of the ag- 
gressive driver. Figure 3 shows the evolution of the density of stopped cars 
in the lanes when starting from the initial homogeneous state. In VDRM , 
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no stopped cars exist in the circuit because all cars respect the safety criteria 
of lane changing. However, in VDRM4, the density increases with time until 
it reaches a stationary value. Stopped cars act as perturbations for the free 
flow region and as such trigger the breakdown of the high flow states. 
In figure 4, we show the cluster size distribution in the steady state of VDRMi 
for different lane-changing probability p c h- The cluster means here a string of 
successive stopped cars in a single lane of the two-lane model, i.e. we are con- 
sidering only compact jams in a single lane. As the symmetric lane-changing 
rules are considered here, the cluster sizes distribution in the two lanes must 
be equal. We observe from figure 4 the bimodal nature of the cluster size 
distribution as p c h — 1. Large clusters appear in the lanes but there are 
by far many more small-sized clusters than large ones. Furthermore, with 
decreasing p c h, the probability of small clusters increases while that of large 
clusters diminishes. If p c h = 0, which corresponds to the single lane VDR 
model, almost all cars are congested in one large cluster with the exception 
of a few isolated cars. This is the well known phase separated state. 
In this section, we shall investigate the relaxation dynamics of VDRMjv a for 
different values of N a and p c h, when starting from the homogeneous initial 
condition. This is done by plotting the time evolutions of the flow and com- 
puting their relaxation times. Hence, the greater is the number of aggressive 
drivers, the faster is the system relaxation (Fig. 5). When decreasing the 
probability of lane changing p ch , one sees that the equilibration is delayed. 
Indeed, in this case, the abrupt lane changes of aggressive cars become less 
frequent and the number of stopped cars becomes small (Fig. 6). 
To study numerically the relaxation time corresponding to an observable A 
we shall use the nonlinear relaxation function [TO] : 

4>{t) = [A(t)-A{o O )]/[A{0)-A{oo)] (1) 
The corresponding nonlinear relaxation time 

POD 

T = <j)(t)dt. (2) 
Jo 

The condition that the system is well equilibrated is 

t Mo > r (3) 

where M is the number of Monte Carlo steps that have to be excluded in 
the averaging of the observable A. In figure 7, we plotted the variation of 
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the relaxation time r of the observable J near the limit p c h — > 0. As a result, 
the relaxation time is found to diverge as p c h —> 0. Moreover, we see that r 
follows a power law behavior of the form, 

r oc Pch ( 4 ) 

Except for some minor fluctuations, the dynamic exponent (3 remains un- 
changed when varying the number of aggressive drivers present in the circuit. 
For example, (3 « 0, 9798 ± 0, 0290 for N a = 1 and (3 « 0, 9801 ± 0, 0315 for 
N a = 3. Assuming that the parameter p c h is rate of transition for the dy- 
namics of the model, and as it was demonstrated in Ref. [UJ, the exponent 
(3 is expected to be theoretically equal to one. 

In summary, we have shown that the presence of aggressive drivers in the 
circuit breakdowns the state of high traffic flow. In theses systems neither 
phase separation nor metastability can occur and a new stationary state takes 
place. Indeed, the abrupt lane changing of aggressive drivers force the suc- 
ceeding cars on the destination lane to decelerate enough; leading therefore 
to the occurrence of stopped cars and then the formation of jams. 
In the NS model, the cluster sizes distribution decreases exponentially while 
in VDRM , it should depend on the initial state. Yet, for some density p in 
the hysteresis region (p x < p < p 2 ), no clusters appear in the metastable ho- 
mogeneous state. However, a big cluster persists in the phase separated state. 
Nevertheless, in VDRMjv„ where N a ^ 0, the stationary state is composed 
by small and big clusters. The distribution of these clusters should depend 
on the lane-changing probability p c h- As p c h decreases the relaxation time 
t of the system increases and diverges at the limit p c h — ► 0. The relaxation 
behaviour follows a power law behavior of the form, r cx , (3 = 1). 
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Figures captions 



Figure 1. Illustration of the variation of the flow J in the two- 
lane VDR traffic models as a function of the density of cars and 
for different values of the discarded time td% s (p c h — 0.10). 
Figure 2. Time evolutions of the flow for homogeneous and 
megajam initial states (p = 0.12 and p c h = 0.10). 
Figure 3. Time evolution of the density of stopped cars when 
starting from the initial homogeneous state (p = 0.12, N a = 1 
and p ch = 0.10). 

Figure 4. The cluster size distribution in the steady state for 

different lane-changing probability p c h (p = 0.12 and N a = 1). 

Figure 5. Time evolutions of the flow and the density of stopped 

cars, when starting from homogeneous initial configuration, for 

several values of N a (p = 0.12 and p ch = 0.10). 

Figure 6. Time evolutions of the flow and the density of stopped 

cars, when starting from homogeneous initial configuration, for 

several values of p c h (p = 0.12 and N a — 1). 

Figure 7. Variations of the relaxation time r near the limit 

Pch^O (p = 0.12). 
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